Inverse Seesaw Neutrino Signatures at LHC and ILC 



We study the collider signature of pseudo-Dirac heavy neutrinos in the inverse seesaw scenario, 
where the heavy neutrinos with mass at the electroweak scale can have sizable mixings with 
the Standard Model neutrinos, while providing the tiny light neutrino masses by the inverse 
seesaw mechanism. Based on a simple, concrete model realizing the inverse seesaw, we fix 
the model parameters so as to reproduce the neutrino oscillation data and to satisfy other 
experimental constraints, assuming two typical fiavor structures of the model and the different 
types of hierarchical light neutrino mass spectra. With the fixed parameters in this way, we 
analyze the heavy neutrino signal at the LHC through tri-lepton final states with large missing 
energy and at the ILC through a single lepton plus di-jet with large missing energy. We find 
that in some cases, the heavy neutrino signal can be observed with a large statistical significance 
via different flavor charged lepton flnal states. Therefore, we can not only discover the heavy 
neutrinos in the future but also obtain a clue to reveal the origin of the small neutrino mass 
and flavor mixing. 
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1 Introduction 



The current experimental results on the neutrino oscillation phenomena [T], including the recent 
measurements of the so-called reactor angle [21 |3l IH O [6], have established the existence of 
neutrino masses and flavor mixings, which require us to extend the Standard Model (SM). The 
seesaw extension of the SM [7] is the most simple idea for explaining the very small neutrino 
masses naturally, where the SM-singlet heavy right-handed Majorana neutrinos induce the 
dimension five operators leading to very small light Majorana neutrino masses (the seesaw 
mechanism) . The seesaw scale varies from the intermediate scale to the electroweak scale as we 
change the neutrino Dirac Yukawa coupling (Yd) from the scale of top quark Yukawa coupling 
(Yd ~ 1) to the scale of electron Yukawa coupling (Yd ~ 10~^). 

In high energy collider experimental point of view, it is interesting if the heavy neutrino 
mass lies at the TeV scale or smaller, because such heavy neutrinos could be produced at 
high energy colliders, such as the currently ongoing Large Hadron Collider (LHC) and the 
International Linear Collider (ILC) being projected as energy frontier physics in the future. 
However, since the heavy neutrinos are singlet under the SM gauge group, they obtain the 
couplings with the weak gauge bosons only through the mixing via the Dirac Yukawa coupling. 
For the seesaw mechanism at the TeV scale or smaller, the Dirac Yukawa coupling is too small 
(Yd ~ 10~^ — 10~^) to produce the observable amount of the heavy neutrinos at the colliders. 

There is another type of seesaw mechanism so-called the inverse seesaw [8] , where the small 
neutrino mass is obtained by tiny lepton-number- violating parameters, rather than the suppres- 
sion by the heavy neutrino mass scale in the ordinary seesaw mechanism. In the inverse seesaw 
scenario, the heavy neutrinos are pseudo-Dirac particles and their Dirac Yukawa couplings with 
the SM lepton doublets and the Higgs doublet can be even order one, while reproducing the 
small light neutrino masses. Thus, the heavy neutrinos in the inverse seesaw scenario can be 
produced at the high energy colliders through the sizable mixing with the SM neutrinos. 

In this paper, we study the inverse seesaw scenario and the heavy neutrino signatures at 
the LHC and ILC. For the concreteness of our inverse seesaw scenario as well as the stabil- 
ity of the electroweak scale, we consider a simple realization of the inverse seesaw mechanism 
in the context of the next-to-minimal supersymmetric SM (NMSSM) proposed in |9]. The 
model parameters are fixed so as to reproduce the neutrino oscillation data as well as other 
experimental constraints such as precision measurements of the weak gauge boson decays and 
lepton-fiavor-violating decays of charged leptons. We consider two typical cases in fitting neu- 
trino oscillation data. One is that the flavor structure among light neutrinos originates from the 
flavor structure of the neutrino Dirac Yukawa couplings. In the other case, the Dirac Yukawa 
couplings are flavor-blind, and the flavor structure among light neutrinos is provided by the 



1 



lepton-number-violating parameters. Assuming the TeV scale mass for sparticles, we concen- 
trate on the production of the heavy neutrinos with mass less than 1 TeV at the LHC and ILC, 
and calculate the number of signal events. The heavy neutrino signals depend on the origin of 
the flavor structure in the model and the types of the hierarchical light neutrino mass spectra. 
We find that in some cases, the signal of the heavy neutrino productions can be observed with 
a large statistical significance for the final states with different charged lepton flavors. In the 
future we can not only discover the heavy neutrinos but also reveal the flavor structure of the 
model which should be compared with the neutrino oscillation data. 

This paper is organized as follows. In Sec. 2, we introduce a model for the inverse seesaw in 
the context of NMSSM. In Sec. 3 we give the explicit formulas of the heavy neutrino production 
cross sections at the LHC and ILC and of the partial decay widths, which are used in our 
numerical analysis. In Sec. 4, we fix the model parameters so as to satisfy the experimental 
constraints, assuming two typical flavor structures of our model and two types of hierarchical 
neutrino mass spectra. With the fixed parameters, the signal of the heavy neutrinos at the 
LHC and ILC are investigated in Sec. 5. Sec. 6 is devoted for conclusions. 

2 Inverse Seesaw in NMSSM 

As a simple realization of the inverse seesaw mechanism, we consider an extension of the 
NMSSM [To]. In the NMSSM, we introduce one gauge singlet chiral superfield S with even Z2 
matter parity through the following superpotential terms: 

W D \SH,,H, + (1) 

where A and k are the dimensionless constants, Hu, are the MSSM Higgs doublets. We 
assume that the hidden sector breaks supersymmetry (SUSY) and induces soft SUSY breaking 
terms for the MSSM scalers and gauginos at the TeV scale, by which the scalar S and the MSSM 
Higgs doublets develop the vacuum expectation values (VEVs). The VEV of S generates the 
MSSM /i term: /i = X{S). 

We extend the NMSSM by introducing new particles and a discrete symmetry [9]. The 
charge assignments for particles relevant to the inverse seesaw mechanism are shown in Table 1. 
Here, and Nj are the MSSM singlet particles, heavy neutrinos of j-th generation, and 
CO = e*^'^/^. There are several possibilities for the Z3 charge assignments, and see [9] for complete 
lists. 

The renormalizable superpotential involving the new particles and being consistent with all 
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Table 1: The charge assignments of the NMSSM superfields. 

the symmetries is given by 

W D YijLiH^m + i^^SNiNj + muN.N^. (2) 

Without loss of generality, we have worked out in the basis where the charged lepton Yukawa 
matrix and m are diagonalized. When the VEVs of S and are developed, we rewrite the 
superpotential as 

W D u'^mnN" + ^N^fiN + N^mN", (3) 

where we have used the matrix notation for generation indeces, v is the MSSM neutrino chiral 
superfield, mn = Yv sin 13 /\/2 with v = 246 GeV is the neutrino Dirac mass matrix, and 
fi = Xn{S). For m larger than the electroweak scale, we integrate out the heavy fields A^^ and 
under the SUSY vacuum conditions, 

dW 

^ N' = -m-^fiN, 



dN 



dW 

—- = 0-^N = -{mom-^fv, (4) 

and we arrive at the effective superpotential at low energies, 

Wcs = ^i^^ [(m£,m^^)/i(m£,m~^)^] u. (5) 

Note that the light Majorana neutrino mass matrix, mj, = {mijm~^)fi{m£,m~^)'^ , is propor- 
tional to /i, so that tiny neutrino masses can be realized by a small fi even for both m and 
rriD being the electroweak scale. This is the inverse seesaw mechanism, where the tiny neutrino 
mass corresponds to the breaking of the lepton number by the tiny fi values. 

Note that the heavy fields being integrated out also have an impact on the Kahler potential. 
Substituting the SUSY vacuum conditions into the canonical Kahler potential for the heavy 
fields, / d^e{N^N + N^^N"), we obtain 

^cff = [(m£,m~^)*(m£,m"^)^] z/ H , (6) 
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where the elhpsis denote higher order terms. Following the electroweak symmetry breaking, 
this dimension six operator induces flavor-dependent corrections to the kinetic terms of the 
left-handed neutrinos [TTj . 

Assuming m/^m"^ ^ 1, we can express the flavor eigenstates (z/) of the light Majorana 
neutrinos in terms of the mass eigenstates of the light (z/m) and heavy {Nm) Majorana neutrinos 
such as 

zz-ATz/^ + W^, (7) 

where 

7^ = mDm-\ Af= (^1- ^e^ Umns (8) 

with e = TZ*TZ^ , and Umns is the usual neutrino mixing matrix by which the mass matrix rrii, 
is diagonalized as 

Umns^uUmns = diag(mi, mg, mg). (9) 

In the presence of e, the mixing matrix A/" is not unitary. Using the mass eigenstates, the 
charged current interaction in the Standard Model is given by 

Ccc = -^W.eYPL [Mv^ + W„) + /i.e., (10) 

where e denotes the three generations of the charged leptons in the vector form, and Pl = 
|(1 — 75) is the projection operator. Similarly, the neutral current interaction is given by 

= -T^Z^ [v;;;YPLi^^^>m + N:::YPLin^'^)N^ + {i7;;i^^PL{M^n)N^ + h.c.]] , (11) 

ZCyj 

where = cosOw is the weak mixing angle. Because of non-unitarity of the matrix A/", 
A/'^A/' 7^ 1 and the flavor-changing neutral current occurs. 



3 Productions and decays of heavy neutrinos at colliders 

In the previous section, we have found the charged and neutral current interactions involving 
the heavy neutrinos. For detailed analysis, we need the information of the mixing matrices, 
A^ and TZ. In the next section, we will flx all the elements of the matrices by considering the 
current experimental results. Before the analysis for flxing the parameters, in this section we 
give the formulas for the production cross sections and the partial decay widths of the heavy 
neutrinos in the limit of one generation and M = TZ = 1. 
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3.1 Production cross section at LHC 

At the LHC, the heavy neutrinos can be produced through the charged current interactions by 
the s-channel exchange of the W boson. The main production process at the parton level is 
ud — e~^Ni (and ud — ?■ e~A'^i) and the differential cross section is found to be 

- (3.89 X io» pb) X JLi+if! f AVj!(^;.-^?i?+^r:?'. d^) 



rfcos^ ' ■ ' 327rs s \2 J \3 J 4 (s - ml^y + m^^Tl, ' 



where is the center-of-mass energy of the colhding partons, M is the mass of A'^i, and 
l3 = {s-M^)/{s + M^). 

The total production cross section at the LHC is given by 



aiHC = [ d\fl I dcos9 I dx—^ — fu{x,Q)fd(—^ ,Q) 

J J Js/El^^s ^^CM5 \xEcMS J 



'^/EhMs ^^CMS \xEcMs J dcosO 

+ {u^u,d^d), (13) 



where we have taken Eqms = 14 TeV for the center-of-mass energy of the LHC. In our numerical 
analysis, we employ CTEQ5M [12] for the parton distribution functions for w-quark (/„) and 
J-quark (/j) with the factorization scale Q = The total cross section as a function of M is 
depicted in Fig. 1. Since we have fixed A/" = 7?. = 1 in this analysis, the resultant cross section 
shown in Fig. 1 corresponds to the maximum values for a fixed M. 

There are three main modes for the heavy neutrino decays: Ni — j- e'W'^, uiZ, uih. The 
corresponding partial decay widths are respectively given by 



r(iVi ^ e-w+) 
T{Ni uiZ) 



g' {M^-ml,f{M^ + 2ml) 
QAtt M^ml, 

(M^ - ml)\M^ + 2m|; 
1287rc2 M%4 



Recently, the ATLAS [13] and the CMS [11] have announced the discovery of a scalar particle 
with mass 125-126 GeV, which is most likely the SM Higgs boson. According to this discovery, 
we use vTLh = 125.3 GeV [T^ in the following analysis. 



3.2 Production cross section at ILC 

The ILC can produce the heavy neutrino in the process e+e~ — ViNi through t and s-channel 
exchanges of W and Z bosons, respectively. The total differential production cross section for 
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Figure 1: The total production cross section of the heavy neutrino as a function of its mass 
at the LHC with ^/s = 14 TeV (sohd hne). As a reference, the production cross section at the 
LHC with ^/s = 8 TeV is also plotted (dashed line). 
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The total production cross section for the process e'^e — t- uiNi for the ILC with 
GeV and 1 TeV, respectively, are shown in Fig. 2. Since we have fixed M = TZ = 



(16) 

s = 500 
1 in this 
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analysis, the resultant cross section shown in Fig. 2 corresponds to the maximum values for a 
fixed M. 
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Figure 2: The total production cross section of the process e+e — ?■ viNi at the ILC with 
a/s = 500 GeV (solid line) and a/s = 1 TeV (dashed line). 

4 Fixing the matrices M and IZ 

The elements of the matrices M and TZ are constrained by the experimental results [IHl [T6l fT7] . 
We begin with the current neutrino oscillation data. Recently non-zero reactor neutrino angle 
6*13 has been observed in several experiments, such as T2K [2], MINOS [3], Double CHOOZ 
Daya Bay [5] and RENO [6], and their results are consistent with each other. Together with 
other oscillation data, all neutrino oscillation parameters except the Dirac CP-phase, two mass 
squared differences and three mixing angles have been measured in some precision. By using the 
data, we fix the neutrino mixing matrix elements. In the following analysis, we use sin^ 26*13 = 
0.092 which is the center value reported by Daya Bay experiments [5]. In their analysis, the other 
oscillation data have been fixed as sin^ 2^12 = 0.87, sin^ 26*23 = 1-0, Am^g = ml— ml = 7.6x10"^ 
eV^, and Am23 = {m"^ — m^l = 2.4 x 10"'^ eV^, and we adopt the same values. Then, the 
numerical values of neutrino mixing matrix elements are explicitly given by 




(17) 
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where we have fixed all the CP-phases to be zero, for simplicity. In fact, we have checked the 
production cross sections of the heavy neutrinos, which will be discussed in the next section, 
weakly depend on the CP-phases, and the introduction of the CP-phases little changes our 
conclusions. 

For the neutrino mass spectrum, we consider both the normal hierarchy (NH) and the 
inverted hierarchy (IH). The lightest mass eigenstate is assumed to be very light and its mass 
is approximated as 0. Thus, in the NH case, the diagonal mass matrix is given by 



Dnh = diag ( 0, -\/ Amfj, \/ Amfg + Am^g 



while in the IH case 

DiH = diag (^y^Amlg- Am?^, ^J Amj^, 0^ . (19) 

In order to make our discussion simple, we assume the degeneracy of the heavy neutrinos in 
mass such as m = Ml with the unit matrix 1, so that the neutrino mass matrix is simplified 
as 

= TZfiTZ^ = -^rriDfim^. (20) 

Moreover we consider two typical cases for the fiavor structure of the model: (i) fi is also 
proportional to the unit matrix, /i — )■ /il. In this case, the fiavor structure of m,^ is provided 
by a non-diagonal rriD. We call this case Flavor Non-Diagonal (FND) case, (ii) The other 
case is what we call Flavor Diagonal (FD) case, where rriD is proportional to the unit matrix, 
rriu — )■ rriul and thus the fiavor structure is encoded in the 3x3 matrix fi. 

In the FND case, we consider two generations of Nj and Nj with j = 1,2, so that 

^•^ = '^^Dm^ = UIjnsDnh/ihUIi^s- (21) 

From this formula, we parameterize the neutrino Dirac mass matrix as 

M ^ 

mo — -j=UmnsV ^NH/iH-i (22) 



where the matrices denoted as ^/D^smiiH are defined as 



^/D^ = I y/Amj^ I , = I a/A^ I . (23) 




Am^g + Aml^J V 
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Note that in the case with two generations of Nj and Nj, the hghtest mass eigenvalue is exactly 
0. On the other hand, in the FD case, we have 

= (^) = UlmsDNH/iHUliNS- (24) 

Next we consider the experimental bound due to the non-unitarity of the mixing matrix M . 
We use the constraints obtained by the combined data from neutrino oscillation experiments, 
precision measurements of weak boson decays, and the lepton-flavor- violating decays of charged 
leptons [16], 

/1.002± 0.005 < 7.2x10"^ < 1.6 x lO^^X 
lA^A/'T] = < 7.2 X 10-5 1.003 ±0.005 < 1.3 x 10^^ . (25) 
\< 1.6x10-2 < 1.3x10-2 1.003 ±0.005/ 

Since MM'^ ^ 1 — e, we have the constraints on e such that 

^0.002 ±0.005 < 7.2 X 10-5 < 1.6 x lO'^X 

< 7.2 X 10-5 0.003 ± 0.005 < 1.3 x 10"^ . (26) 

< 1.6x10-2 < 1.3x10-2 0.003 ±0.005/ 

The most stringent bound is given by the (1,2) element which is from the constraint on the 
lepton-fiavor-violating muon decay /i — ?■ exi- For the FND case, we describe e of the form, 

e = j-p'^D^]:) = -UmnsDnh/ihUm^s, (27) 

and we fix the minimum fi value (fimin) so as to give ei2 = 7.2 x 10"^ with the oscillation 
data in Eqs. ([I7D, ([HD and ([19]). We find fimin = 109 eV and 68.6 eV for the NH and IH 
cases, respectively. Here we have used the fact that all parameters are real according to our 
assumption. In this way, we can completely determine the mixing matrix TZ and A/" from Eq. ( 122]) 
by taking /i = /imiru which optimizes the production cross sections of the heavy neutrinos at the 
LHC and ILC. For the FD case, e = {mD/M fl and we simply take e = diag(0.007, 0.007, 0.007) 
as a conservative treatment. 



5 Collider signatures of heavy neutrinos 

Let us now analyze the collider signatures of the heavy neutrinos with the information of TZ 
and A/" determined by the previous sections. In Sec. 3, we have already given the formulas 

^ Recently it has been pointed out [HI [19] that in the SUSY inverse seesaw model, sparticle Z-penguin 
contributions can dominate the lepton-flavor-violating processes, independently of sparticle mass spectrum. 
According to the analysis in Ref. |19| , we have found that the constraint from fi ~ e conversion process is more 
severe than the one from the /i — >■ 67 process for M > 335 GeV. Since we will focus on M = 100 — 150 GeV in 
the following analysis, we use the value of fimin determined from the muon decay constraints. 
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Figure 3: Signal cross sections proving the tri-lepton final states for the FND (left panel) and 
FD (right panel) cases, at the LHC with ^/s = 14 TeV. 

used in our analysis in the limit of 7^ = A/" = 1. It is easy to generalize the formulas with the 
concrete TZ and A/". The production cross section of the i-th generation heavy neutrino at the 
LHC, through the process qq' — )■ iNi {ud — )■ f^A^i and ud — > i^^i) is given by 

a{qq' ^ i^N.,) = aLHc\n^i?, (28) 

where ainc is the cross section given in Eq. f[T^ . Similarly, the production cross section at 
the ILC is 

a(e+e- ^ V^Ni) = aiLcl^J^ (29) 

where ajLc is given in Eq. ( HSj) . and we have used the approximation Af^R ^ Umns^ because 
Ica/?! <^ 1 as discussed in the previous section. The partial decay widths for the process 
Ni — )■ l~W^ /vaZ /vah are obtained by multiplying Eq. f lT^ and the factor |7?.ojp together. 

5.1 Heavy neutrino signal at LHC 

As has been studied in Ref. [20] (see also ||21j for the studies on the left-right symmetric model), 
the most promising signal of the heavy neutrino productions at the LHC is obtained by the 
final state with three charged leptons (£^£^£^ with the total charge ±1) through the process 
qq' — )■ Nl^ followed by — > i^W^ and — i^u. In this work, detailed studies have been 
performed for the signal of the heavy neutrino with a 100 GeV mass, which couples with either 
the electron or the muon. The events were pre-selected for two like-sign charged leptons (ee 
or fifi) to have transverse momentum px > 30 GeV. The decay mode, uZ, followed by 

Z — £~^i~ is rejected by a cut for the invariant mass of the charge neutral di-lepton. After 
elaborate selections, it has been concluded [20] that the heavy neutrino coupling to the muon 
could be observed at the LHC through the tri-lepton final states. 
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ee 




FND (NH) 


0.30 


1.80 


FND (IH) 


8.24 


3.78 


FD 


104.5 


100.2 


SM background 


116.4 


45.6 



Table 2: Number of events at the LHC with y/s = 14 TeV and 30 fb ^ luminosity, for the 
heavy neutrino mass M = 100 GeV. 

In our analysis, we follow the procedure in Ref. [20]. Since we are considering the general 
case with TZ and A/" consistent with the experimental data, the production cross sections of 
the heavy neutrinos are different. Fig. 3 shows the signal cross section providing tri-lepton 
final states with ee or fifi for the FND (left) and FD (right) cases, as a function of the heavy 
neutrino mass. In the left panel, the dashed and solid lines correspond to the NH and IH cases, 
respectively. The upper solid (dashed) line shows the cross sections with ee (/i/i). 

We adopt the same efficiency for the signal events and the SM background events which was 
found in Ref. [20]. The number of events for tri-lepton final states with ee and fifi, respectively, 
are listed on Table 2, for the luminosity 30 fb~^. Unfortunately, the number of events for 
the FND case is found to be too small. This is because the component of TZ is very severely 
constraint to be small by the current experiments. On the other hand, the FD case results 
the large number of signal events. When we naively estimate the significance by S/ \/B, the 
luminosity of 8.0 fb~^ (3.4 fb~^) is enough to achieve 5-cr significance for the ee (/i/i) final states. 

5.2 Heavy neutrino signal at ILC 

The signature of heavy neutrinos at the ILC has been studied in detail based on the realis- 
tic Monte Carlo simulations in Ref. [22]. In the studies, a five-dimensional model with bulk 
right-handed neutrinos [23] is considered and its 4-dimensional effective theory provides the 
Kaluza-Klein tower of the heavy neutrinos having sizable coupling to the weak gauge bosons 
through mixings with the SM light neutrinos. This structure of the couplings between the 
heavy neutrinos and the SM particles in the five-dimensional model is similar to the one in 
our inverse seesaw model. Thus, we apply the results in Ref. [22], in particular, the signal and 
background selection procedure to our model. 

According to Ref. [22], we focus on the two-jets and one isolated lepton signal with large 
missing energy at the ILC: e"''e~ — uN, followed by the decays — £W and W — )■ qq', 
through which the heavy neutrino production cross sections and the heavy neutrino mass can 
be reconstructed. The production cross sections for the process e+e" — )■ uN, followed by the 
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Figure 4: The production cross sections for the process e^e~ — )■ vN , followed by the decays 
— {i = e, /i, r) and W — )■ qq', as a function of the heavy neutrino mass. The upper-left 
panel shows the results for the FND case with ^/s = 500 GeV. The upper-right panel is the 
same as the upper-left panel but for the case with y/s = 1 TeV. The results for the FD case are 
shown in the lower panel for ^/s = 500 GeV (solid) and ^/s = 1 TeV (dashed), respectively. 



decays — )■ £W {£ = e, /i, r) and W — )■ qq', as a function of the heavy neutrino mass are 
depicted in Fig. 4. The upper-left panel shows the results for the FND case with ^/s = 500 
GeV. The dashed lines from top to bottom denote the signal cross sections for i = t, fi 
and e, respectively, in the NH case. The two lines corresponding to i = t and fi are well- 
overlapping. The solid lines from top to bottom denote the signal cross sections for i = e, t 
and /X, respectively, in the IH case. The two dashed lines corresponding to i = fi and r are 
well-overlapping. The upper-right panel is the same as the upper-left panel but for the case 
with ^/s = 1 TeV. The results for the FD case are shown in the lower panel for ^/s = 500 GeV 
(dashed) and ^/s = 1 TeV (solid), respectively. Here i is either e, fi or r. 

For the ILC with ^/s = 500 GeV and the luminosity C = 500 fb~^, the signal and background 
events are listed on Table 2, Here, the final state of one electron and two jets with missing energy 
from anti- neutrinos is considered, and we have adopted the efficiencies found in Ref. [22], 
for M = 150 GeV. The main backgrounds are euW — ?■ euqq and WW — ?■ iuqq, which are 
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dramatically reduced by the selection using an isolated-electron track with a requirement of its 
energy range, 10 GeV< E^. < 200 GeV, the requirement of the reconstructed di-jet mass to be 
consistent with W hypothesis etc (see [22] for details). The signal and background events for 
y/s = 1 TeV and the same luminosity £ = 500 fb~^ are listed on Table 3. For completeness, 
we have also listed the signal events (without cuts) for the case i = fi and r. 





Events before cuts 


Events after cuts 


FND (NH) 


594 


403 


END (IH) 


11504 


6542 


FD 


58406 


39537 


SM background 


3210500 


23346 


FND (NH, i = n) 


4067 




FND (NH, e = T) 


4257 




FND (IH, i = ij) 


6050 




FND (IH, i = T) 


6077 





Table 3: The number of events at the ILC with ^/s = 500 GeV and the luminosity 500 fb~^, 
for the heavy neutrinos with mass 150 GeV. We have adopted the efficiencies for the signal and 
the SM background found by the realistic Monte Carlo simulations in [22] . 





Events before cuts 


Events after cuts 


FND (NH) 


776 


248 


FND (IH) 


15038 


3725 


FD 


76346 


24384 


SM background 


5476408 


10500 


FND (NH, i = fi) 


5316 




FND (NH, i = r) 


5565 




FND (IH, i = fi) 


7909 




FND (IH, i = T) 


7944 





Table 4: The same as Table 3, but ^/s = 1 TeV. 

In both y/s = 500 GeV and 1 TeV, the signal to background ratio is quite large for the 
FND (IH) and FD cases, while the significance is around 2.5-cr for the FND (NH) case. We can 
expect a similar significance for the i = fi case, and the heavy neutrinos can be observed with a 
large significance for all the FND and FD cases. In Ref. [22], the i = t case is also analyzed in 
detail. In this case, the signal — ?■ veW{W — )■ qq') is considered as the background, and the 
analysis depends on the number of the signal events and hence, we cannot simply adopt the 
results in Ref. [22]. However, since the main backgrounds are evW — )■ evqq and WW — )■ Ivqq' 
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also for this case, we can expect that the efficiency for our case is similar to the one obtained in 
Ref. |22], which is roughly the same as in the I = e case. Thus, we expect a large significance 
for the signal events for all FND and FD cases. 

6 Conclusions 

We have studied the inverse seesaw scenario and the signature of the pseudo-Dirac heavy 
neutrino production at the LHC and ILC. In the inverse seesaw scenario, the light neutrino 
masses are realized by small lepton-number-violating parameters and hence the SM singlet 
neutrinos have sizable Dirac Yukawa couplings with the SM lepton doublets and Higgs doublet 
even for their mass scale being at the TeV scale or smaller. As a result, the heavy neutrinos can 
be produced at the LHC and ILC. Based on a concrete model realizing the inverse seesaw in the 
context of the NMSSM, we have fixed the model parameters so as to satisfy the experimental 
results such as the neutrino oscillation data, the precision measurements of the weak boson 
decays, and the lepton-flavor-violating decays of charged leptons. We have considered two 
typical cases for the neutrino flavor structures of the model, namely, the FND and FD cases. 
With the fixed parameters, we have calculated the production cross sections of the heavy 
neutrinos at the LHC and ILC. 

For the LHC with y/s = 14 TeV, we have analyzed the productions of the heavy neutrinos 
with a degenerate 100 GeV mass, providing the tri-lepton final states with the like-sign electrons 
or muons. After imposing suitable cuts, we have found that the 5-cr statistical significance of 
the signal events over the SM background can be achieved even for the luminosity less than 10 
fb~^ in the FD case. On the other hand, the production cross sections in the FND case is too 
small to observe the heavy neutrino signal. 

We have also studied the heavy neutrino production at the ILC with ^/s =500 GeV-1 TeV, 
where the final state with a single, isolated electron, and di-jet and large missing energy is 
considered. For the luminosity a/s = 500 fb~^, we can obtain clear signatures of the heavy 
neutrinos with mass 150 GeV for the IH mass spectrum in the FND case and the FD case. On 
the other hand, we have obtained around 2.5-0" significance for the NH mass spectrum in the 
FND case. Since we can expect the similar efficiencies of the signal and SM background for 
the final states with different lepton flavors, muon or tau, the heavy neutrinos can be detected 
with a large statistical significance in the modes for all FND and FD cases. 

The flavor structure of the model reflects the flavor dependent signal events from the heavy 
neutrino productions. Once the heavy neutrinos are observed in the future, the flavor depen- 
dence of the signal events gives us a clue to reveal the origin of the neutrino masses and flavor 
mixings. Note that observations of the SM singlet heavy neutrino signal already imply that a 
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mechanism of the neutrino mass generation is not due to the conventional seesaw mechanism. 

Finally we comment on the current bound of the heavy neutrino production at the LHC 
The ATLAS experiment [23] has reported their results on the search for the heavy neutrinos 
based on the production through effective four-fermion operators |25]. The vector operator of 
{d'~f^u){N'~ffj,i)/A'^ is relevant to our case. The final states with iijj {£ = e or /i) have been 
analyzed as a signal of the heavy neutrino production, followed by the decay — iW, W — )■ jj. 
From the data corresponding to an integrated luminosity of 2.1 fb~^ at y/s = 7 TeV, the ATLAS 
experiment has set the lower bound on the cutoff scale A as a function of the heavy neutrino 
mass > 200 GeV. For example, it is found that A > 2.8 TeV for M = 200 GeV. We interpret 
this result to the upper bound on the heavy neutrino production cross section through the 
four-fermion operator as cr{qq' -^■ ^A^) < 24.0 fb. In the FD case, we find cr{qq' — > iN) ^ 3.77 
fb and therefore, the parameter region we have examined in this paper is consistent with the 
current LHC results. 
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